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l> ■ ABSTRACT 

It is generally thought that most of the spin-down power of a pulsar is carried away 
in an MHD wind dominated by Poynting flux. In the case of an oblique rotator, a signif- 
icant part of this energy can be considered to be in a low-frequency wave, consisting of 
stripes of toroidal magnetic field of alternating polarity, propagating in a region around 



the equatorial plane. Magnetic reconnection in such a structure has been proposed as a 
mechanism for transforming the Poynting flux into particle energy in the pulsar wind. 
We have re-examined this process and conclude that the wind accelerates significantly 
in the course of reconnection. This dilates the timescale over which the reconnection 
^ process operates, so that the wind requires a much larger distance than was previously 

thought in order to convert the Poynting flux to particle flux. In the case of the Crab, 
the wind is still Poynting-dominated at the radius at which a standing shock is inferred 
from observation. An estimate of the radius of the termination shock for other pulsars 
implies that all except the milli-second pulsars have Poynting-flux dominated winds all 
the way out to the shock front. 

Subject headings: pulsars: general — pulsars: individual (Crab) — MHD — stars: winds 
and outflows — plasmas — waves 



1. Introduction 

The diffuse synchrotron radiation from the Crab Nebula has now been observed in great detail 
in several wavebands, e.g., Hester et al. (1995). Although by far the best observed example, the 



Crab is just one of several pulsars surrounded by such a nebula (Arons 1996), whose ultimate 
source of energy is almost certainly the rotational kinetic energy of the central neutron star (Pacini 
1967). However, despite intensive efforts, it is still not known how this energy is transferred from 
the neutron star to the radiating, relativistic electrons. 

It is widely accepted that pulsars emit an electron-positron plasma, which carries away energy 
in the form of an ultra-relativistic magnetized wind, together with large amplitude waves (Rees 
& Gunn 1974). At a shock front, located, in the case of the Crab Nebula, some 10 17 cm from 
the neutron star, this energy is released into the relativistic electrons responsible for the observed 
radiation. The most serious problem with this picture is that close to the strongly magnetized 
neutron star the energy must be carried mostly by electromagnetic fields as Poynting flux (Michel 
1982). But, in order to produce the radiating electrons, the energy flux at the shock front must 
be carried mainly by the particles (Rees & Gunn 1974; Kennel & Coroniti 1984; Emmering & 
Chevalier 1987). As has been pointed out by several authors -- and emphasized recently by 
Chiueh, Li & Begelman (1998) and Bogovalov & Tsinganos (1999) — in an ideal, ultra-relativistic 
MHD wind, there is no plausible way of converting the Poynting flux into particle energy flux. 
Several effects have been investigated in attempts to overcome this difficulty, ranging from rapid 
expansion in a magnetic nozzle (Chiueh, Li & Begelman 1998) to non-ideal MHD effects in a 
two-fluid (electron and positron) plasma (Mestel & Shibata 1994; Melatos & Melrose 1996). It 
has also been suggested that a global non-axisymmetric instability of the nebular plasma alleviates 
the problem by enabling the conversion of magnetic to particle energy (Lyubarsky 1992; Begelman 
1998). However, perhaps the most promising explanation to date has been that of reconnection in a 
striped pulsar wind, as advanced by Michel (1982), Coroniti (1990) and Michel (1994). Here, we 
reanalyze this process. Both Coroniti and Michel effectively assumed the wind maintained constant 
speed during reconnection. We find that the hot plasma in the current sheets where reconnection 
takes place performs work on the wind, accelerating it substantially. As a result, the reconnection 
rate is much slower than hitherto claimed. In the case of the Crab Nebula, we conclude that 
reconnection is not capable of converting a significant fraction of the energy flux before the wind 
reaches the position at which the shock front is encountered. 

The conversion of Poynting flux to kinetic energy is also a central problem in a certain class 
of models of gamma-ray bursts. It has been proposed that the emission of a large amplitude elec- 
tromagnetic wave by a black hole or neutron star(s) underlies the burst phenomenon, that the 
gamma-rays are generated by dissipative processes in the wave, and that the X-ray and lower fre- 
quency emission is produced at an outer shock front (Usov 1992, 1994; Blackman & Yi 1998). The 
calculations we present can be rescaled to this situation and may offer a mechanism for accelerating 
the wind to the very high Lorentz factors required before dissipation sets in. 

The paper is organized as follows: in Section 2 we describe the problem as formulated by 
Coroniti (1990) and Michel (1994), present estimates of the effect of reconnection and explain 
in physical terms the reasons for our new result. The perturbation method we employ is a short 
wavelength approximation, which uses as a small parameter the ratio of the radius of the light 



cylinder to the actual radius. This is presented in outline in Section 3, relegating most of the 
algebra to the Appendix, but pointing out the differences between our method and that of Coroniti 
(1990). An analytic asymptotic solution to the system is given in Section 4, followed by results of a 
numerical integration of the equations and a discussion of the limitations of our approach. Finally, 
in Section 5, we summarize our conclusions and their implications, especially for our understanding 
of the Crab Nebula. 



2. The striped wind 

Despite the presence of a plasma, the total power lost by a rotating magnetized neutron star 
may be estimated using the formula for a magnetic dipole rotating in vacuum Gunn & Ostriker 
(1969), Michel (1982). However, the presence of plasma changes the physical picture significantly: 
for example, even an axisymmetric rotator surrounded by plasma loses energy by driving a plasma 
wind, in contrast to an aligned magnetic dipole in vacuum. Furthermore, although a strong vacuum 
electromagnetic wave readily loses energy to particles (Ostriker & Gunn 1969; Asseo, Kennel & 
Pellat 1978; Melatos & Melrose 1996), the presence of plasma is likely to prevent the formation of 
such waves and restrict dissipation to shock fronts or regions in which the magnetic field undergoes 
reconnection. 

In the case of an oblique rotator, the energy lost in the wind can be regarded as shared 
between an axisymmetric component of the Poynting flux and one due to MHD waves, the ratio 
being determined by the angle between the magnetic and rotational axes. Michel (1971) pointed 
out that such waves, which have a phase speed less than that of light, should evolve into regions of 
cold magnetically dominated plasma separated by very narrow, hot, current sheets. The formation 
of this pattern may be imagined as follows: in the axisymmetric case, a current sheet separates 
the two hemispheres with opposite polarity beyond the light cylinder. As the obliquity increases, 
this sheet begins to oscillate about the equatorial plane because the field line at a given radius 
alternates in direction with the frequency of rotation, being connected to a different magnetic pole 
every half-period. Such a picture is observed in Solar Wind. In a quasi-radial flow, the amplitude 
of these oscillations grows linearly with radius and at large distances one can locally imagine quasi- 
spherical current sheets following each another and separating the stripes of magnetized plasma 
with opposite polarity. Coroniti (1990) called this picture a striped wind. Recently Bogovalov 
(1999) has found an exact solution for the oblique split monopole case which has precisely this 
structure. 

It was noticed by Usov (1975) and Michel (1982) that these waves must decay at large 
distances, since the current required to sustain them falls off as r _1 . This is slower than the fall-off 
in the available number of charge carriers, which goes as r~ 2 . Coroniti (1990) considered the 
reconnection process in a striped wind and came to the same conclusion. In the case of a highly 
oblique rotator, both Coroniti and Michel agreed that the MHD wind of the Crab pulsar could be 
transformed by this mechanism from one dominated by Poynting flux to one dominated by particle 



kinetic energy flux well within the radius of the termination shock. 

The pulsar wind beyond the light cylinder is quasi-spherical (Chiueh, Li & Begelman 1998; 
Bogovalov & Tsinganos 1999), and, as the magnetic field is predominantly toroidal, it scales roughly 
as 

B = £ L -, (1) 

r 

where tl = c/Q, = cP/2ir is the light-cylinder radius, O the angular velocity of the neutron star P 
its period of rotation and B^ the magnetic field strength at j-l- Near the light cylinder, the poloidal 
and toroidal components are comparable, but within it the poloidal field is nearly dipolar, so that 
one can estimate 

*L = 4 » SSC (2) 
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where [i = 10 30 /U3o G-cm 3 is the magnetic moment of the star and P is given in seconds. The density 
n in the quasi-spherical flow, measured in the rest frame of the star, decreases approximately as 
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n = nL W) ' (3) 



since in a steady, spherical, relativistic flow r 2 nv is constant and the speed v is close to c. The 
density is conveniently normalized by the Goldreich- Julian charge density pgj : 



u-l = Kp G j/e 
2irec 



(4) 



where k is the multiplicity coefficient. This quantity is rather uncertain but generally expected to 
be large: k ~ 10 3 - 10 4 (Arons 1983). 



The Poynting flux may be estimated as 



c J%r±V 



w = 17 if ; (5) 

and the ratio of the Poynting flux to the kinetic energy flux, called the magnetization parameter 
a, is 

W 

(6) 
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At the light cylinder this quantity takes on the value 
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where wl = eB^/mc is the gyrofrequency at the light cylinder. The ratio uj-l/Q is large for all 
pulsars. For example, in the case of the Crab, where ^30 ~ 5, we have wl/Q ~ 10 11 . 

The speed of a fast magnetosonic wave propagating perpendicular to the magnetic field in a 
magnetically dominated plasma corresponds to a Lorentz factor 7f ms = \fa [e.g., Kirk & Duffy 
(1999)]. In pulsar models, plasma is ejected at Lorentz factors of 10 2 to 10 3 , which is higher than 
7mfs f° r most pulsars. In the case of the Crab, however, these values correspond to a trans-alfvenic 
speed (7L ~ \f&V) at the light cylinder. At some point beyond the light cylinder, but before the 
onset of dissipation, we expect the wind to establish a supersonic ideal MHD flow, in which the 
value of a and 7 are constant. In the following, these values are indicated by the subscript 'L', even 
though, strictly speaking, they may not be achieved at the light cylinder itself. Using the above 
estimates, we then find that, for the Crab, k = 10 4 corresponds to ox ~ 3 x 10 4 and 7l ~ 200. At 
the light cylinder, almost all the energy is carried by the electromagnetic field. If this energy were 
completely transferred from Poynting flux into plasma kinetic energy, the Lorentz factor would 
attain the value 



7max = 7LO"L 



(9) 



However, a cold, radial MHD wind does not accelerate, because the outwardly directed pressure 
gradient arising from the gradient of the magnetic field is exactly balanced by the inwardly directed 
tension force exerted by the curved toroidal field lines. In the absence of dissipation, the energy 
flux remains locked in the field as Poynting flux. 

The idealized radial structure of the wind and embedded current sheets corresponding to the 
field of a perpendicular split monopole is shown in Fig. 1. In the proper plasma frame the current 
density I' (current per unit length in the toroidal direction) is, according to Ampere's law, 

Here and in the following we use primed values to indicate quantities in the proper plasma frame, 
e.g., B' = -B/7, n' = n/7 etc. Usov (1975) and Michel (1982) noticed that the current in Eq. (10) 
cannot be maintained to arbitrarily large radius: since the proper current density cannot exceed 
en'c and the width of the sheet is less than half a wavelength, nr^, one can easily see from Eqs.(l), 
(3) and (4) that the striped wind cannot exist beyond the radius 

»"max = 7TK77-L (11) 

[see also (Melatos & Melrose 1996)]. When the velocity of the current carriers approaches the speed 
of light, an anomalous resistance arises and the alternating magnetic field dissipates by heating the 
plasma. Coroniti (1990) considered the dissipation as reconnection through the current sheet, 
whose minimum thickness he took equal to the Larmor radius which a thermal particle of the hot 
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Fig. 1. — Idealized picture of a pulsar wind driven by a rotating, perpendicular, split-monopole 
magnetic field, showing the magnetic field, B (solid line), particle density n (dashed line) and 
plasma pressure p (dotted line) as functions of the phase variable <j> = r/rx,, with origin at an 
arbitrary position in the outer part of the wind (r 3> n,)- The magnetic field and density are 
constant between the current sheets, which are shown shaded. Within the sheets, we assume zero 
magnetic field and constant density. 



plasma in the sheet would have, if it entered the cold magnetized part of the wind. In terms of the 
fraction A of a wavelength 2irr^ occupied by the two current sheets, this condition reads 

T' 

where T" is the characteristic "temperature" of the particles in the sheet, in energy units. Taking 
into account that the sheet is in pressure equilibrium, 

n'T' = — (13) 

one can easily see that the current density j' < en'c, so that Coroniti's condition is, to within a 
factor of the order of unity, the same as that of Michel (1994). 

If, initially, a fraction a of the plasma is concentrated in the sheets, dissipation begins at a 
radius r\ = ax max - For r < r\ the sheet width exceeds the limiting value of Eq. (12), the current 
carriers remain non-relativistic, the conductivity is very high and reconnection does not proceed. 
For r > n, however, the current carriers become relativistic, an anomalous resistivity arises, and 
reconnection brings energy and plasma into the sheet until, eventually, all the magnetic flux is 
destroyed at the radius r = r Tj 



max ■ 



Both Michel and Coroniti estimated the radius r max by substituting into Eq.(ll) the initial 
Lorentz factor tl. The corresponding value is large compared to the light cylinder radius. For the 
Crab, it is of the order of W 7 r^, which is, nevertheless, still well within the radius of the standing 
shock r s ~ 10 9 rL- However, as we show below, the flow accelerates as reconnection proceeds. 
Basically, this is because the hot plasma continues to exert an outwardly directed pressure gradient, 
but the compensating inwardly directed tension force is absent. Equivalently, it is clear that the 
hot plasma performs work during the radial expansion, and this appears as an acceleration of the 
wind. In an accelerating wind, most of the energy is released when the Lorentz factor of the flow 
is roughly at its maximum value, given by Eq. (9). The corresponding radius is 

'"max — T^Tmax^L 
TTU1-L 



r L (14) 
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In the case of the Crab, r max ~ lO 11 ^, which significantly exceeds the radius of the standing shock, 
so that only a small fraction of the wave energy is converted into particle energy before the plasma 
arrives at this shock front. 



3. Equations 

3.1. Local structure of the wave 

The striped wind described above can be regarded as an MHD wind containing an entropy 
wave which moves together with the plasma. At large distances the wave consists locally of spherical 
current sheets separating cold, magnetized stripes of plasma with opposing polarity (see Fig. 1). 
To find the evolution of this wave, we use a two-timescale perturbation approach, assuming that 
the timescale on which it evolves i.e., the timescale on which the reconnection sheet grows, is much 
longer than the period P of the wave. Details of this method are presented in the Appendix. Here 
we restrict ourselves to the most important steps in the argument. 

In the cold magnetized part of the wind the magnetic field B, which is toroidal, is constant on 
the fast timescale and the plasma pressure vanishes. The plasma density n' c (n c in the lab. frame) 
is also constant on the fast timescale. In the hot sheet, the magnetic field vanishes, the pressure 
p' is constant on the fast timescale, as is the plasma density n' h , which, in contrast to previous 
treatments, is not constrained to equal n' c . In the plasma frame the entire pattern is at rest, but 
the plasma speed v, and the corresponding Lorentz factor 7 change on the slow timescale as the 
wave evolves. In the cold part of the wind, the magnetic field is frozen into the plasma, so that in 
spherical polar coordinates we can write for the slow evolution in radius 

d /B\ 

— — = 0. 15 

dr \rn c J 

[see Eq. (A33)]. Although the method is strictly valid only for r 3> rx, we normalize the quantities 
in the wind to their value at r = n,. Using Eq. (4) and interpreting the multiplicity factor as 
referring to the density of electron/positron pairs outside the current sheets, this yields 

4 = -■ <*> 

rn c k 

Following Coroniti (1990), we assume that reconnection keeps the sheet width equal to the limiting 
value given in Eq. (12) and (A32). The condition that the cold and hot parts of the flow are in 
pressure equilibrium (13) leads, together with Eq. (16), to the relation 

A = -^— (17) 



where a useful variable 



Z = 4 = ^ (18) 



has been introduced. In his treatment, Coroniti (1990) set n' h = n' c to write his equation (16) 
the counterpart of our Eq. (17). 



The remaining equations needed to describe the wave evolution are those of conservation of 
particle number and energy (equivalent to the equation of motion in the relativistic formulation) 
and the entropy equation. To zeroth order in r-^/r these confirm that the plasma speed equals the 
pattern speed, the hot and cold parts are in pressure equilibrium and the configuration shown in 
Fig. 1 and described above is stationary on the fast timescale. 



3.2. The continuity equation 

The slow evolution of the system is obtained by averaging over a wave-period, which we denote 
by angle brackets: (•••)■ Conservation of particle number (the continuity equation) gives 

(nvr ) = constant (19) 

This equation is exact. In the short wavelength approximation, one can replace the time average by 
a spatial one, because all slowly varying parameters are constants on the scale of one wavelength. 
This immediately yields 



vr n c 



§ + d-A) 



= C\ = constant (20) 



[see Eq. (A34), noting that in this section we use 7 and v to refer to the zeroth order quantities, 
denoted in the Appendix by 70 and vq\. 



3.3. The energy equation 

After time averaging, the energy equation is 

2 / / 2 EB \ 

r ( w 7 v + — — c ) = constant (21) 

Here w' is the enthalpy density: outside of the sheets w' = mn' c c 2 , whereas within them w' = 
<±p' + mn' h c 2 (= B' 2 /2n + mn' h c 2 ). One can substitute for the electric field using E = vB because, 
by assumption, E and B are nonzero only outside the sheets, where magnetic field is frozen into 
the plasma. 

Replacing once again in Eq. (21) the time average by a spatial one yields 

?/2 ™„/^2\ / r/2 s 



2 2 

7 vr 



B' 2 mn' c c 2 \ . „ A . / , 2 B' 2 \ 



= Co = constant (22) 



[see Eq. (A35)]. 



Coroniti (1990) did not use equation (21) in his analysis. Instead, he assumed the plasma 
both inside and outside of the sheet is strictly stationary in the wave frame so that the sheet edge 
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moved through a constant density plasma. This led him to set n' h = n' c in Eq. (17) (our numbering), 
thus reducing the number of unknowns. He was then able to find an expression for the speed of 
the sheet edge. However, such a picture violates energy conservation, as is apparent from Eq. (21): 
the condition of pressure equilibrium (p 1 = B 2 /8tt^ 2 ) means that a decrease in the enthalpy flux in 
the magnetized part of the flow (= r 2 vB 2 /Air) cannot be balanced by an increase in the enthalpy 
flux in the sheet (= r 2 w'^ 2 v = 4r 2 p'^/ 2 v ) unless there is a velocity (and density) jump across the 
sheet edge. 



3.4. The entropy equation 

The entropy equation requires more care. Following Eq. (21) of Coroniti (1990) and setting 
the ratio of specific heats to 4/3, we write for the full nonlinear equation: 

dp' , p' dn' / _ v \ . . 

3-=- -1—,-ir = (E--B)j (23) 

dt ri dt \ c J y ' 

(see Eq. A4) where the convective derivative d/dt = d/dt + vd/dr. The right-hand side of this 
equation, when multiplied by 7, is the rate of entropy generation in the rest-frame of the plasma, 
which moves with speed v. In the limit of a sharp transition between the magnetized and unmagne- 
tized parts of the flow, the entropy generation term is the product of two singular functions: one for 
the current and another for the electric field in the plasma frame. To find the slow evolution of the 
wave, it is essential to perform the averaging process before inserting this specific representation, 
which is not well-defined at the point at which entropy is generated. In moving from his Eq. (21) 
to (22), Coroniti overlooked this point. As a consequence of this, and of the incorrect expression 
for the expansion speed of the sheets, he came to the erroneous conclusion [in his equation (26)] 
that the wave Lorentz factor remains almost constant during reconnection. 

Following the averaging procedure described in the Appendix (Eqs. A27 to A30), we find [see 
Eq. (A37)]: 

%— (rVA) + 3 1V — (Ap') + 

-g- r b 2 ( 1 - A )p]+- KA - r L Q- r (rv) = 0. (24) 

This equation may be integrated, using the equations of continuity (20) and flux freezing (16): 

(1 - 3Z)n' c r 3 = C 3 = constant (25) 

[See Eqs. (A38-A42)]. 

The system thus consists of five algebraic equations (16), (17), (20), (22) and (25) for the five 
slowly varying unknown functions of radius: n' c , Z, B', A and 7. 



11 



Results 



4.1. Asymptotic solution 



A general solution to this system is difficult to find, and the most straightforward way to 
generate a numerical solution is to revert to integrating the differential forms of the continuity, 
energy and entropy equations. However, it is possible to extract analytically an asymptotic solution, 
valid for A, 7 -2 , c^ <C 1, which is just the regime we are interested in. 

As can be checked a posteriori, the quantity n' c r 3 is an increasing function of r, so that, 
according to Eq. (25), Z — ► 1/3 as r — ► oo. Defining, in accordance with Eq. (7), 

■S /2 /4vr" 



C"L 



mn'c 2 



(26) 



r=r L 



(27) 



we can use the continuity equation (20) to rewrite the energy equation (22) as 

B' 2 

7 2 t>r 2 (l + A)- \--fmc 2 C\ = 7 L mc 2 Ci(l + <r L ) . 

Here we have assumed the initial thickness of the current sheet is vanishingly small. This is 
reasonable, since before the onset of reconnection, the plasma in the sheet expands adiabatically 
and cools, causing the sheet width to decrease. To lowest order in the small parameters, this 
equation (27), together with Eq. (16), and the continuity equation (20) yield the same result: 



jr n' c 



C\ 



(28) 



so that the system is nearly degenerate. Expanding the equations to next order and eliminating 
the leading term gives 
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/L 2 7 V 7L 

At large radius Z — ► 1/3 and for a super- Alfvenic flow (7 3> 7l > s/v), one finds 7/7L 
leading to the asymptotic solution 
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mc 2 
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187TWLrL 
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1/2 



180r 
7ru; L r L 
. -3 

r 



-1/2 



27to;l 

n 



(29) 
3ctlA, 

(30) 
(31) 

(32) 
(33) 



which agrees with our estimate that the maximum value of 7 given in Eq. (9) is attained at the 
radius r max of Eq. (14), where A ~ 1. Note that this solution is independent of the actual value of 
7l, which enters only as a scaling factor for 7. 
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4.2. Numerical solution 

Three physical parameters determine the cold radial MHD wind of a pulsar in the absence of 
reconnection. They are the values at the light cylinder of the magnetization parameter <tl, Lorentz 
factor 7l and the ratio of the particle gyro frequency to the rotation period oj-^/Q.. These are 
related to the multiplicity parameter k by Eq. (7). In the presence of reconnection, the same three 
parameters also uniquely specify the asymptotic solution at large radius Eqs. (30-33). However, 
the full solution requires one additional initial condition, which is the fraction a of plasma which is 
initially present in the current sheets. This quantity determines the radius at which reconnection 
starts, which is, formally, the position at which we impose the initial condition A = 0. The 
perturbation analysis presented above is valid when this radius is large compared to n,. 

To investigate the dependence of the solutions on this initial condition, we have solved the 
system by integrating numerically the three equations continuity, energy and entropy written in 
differential form: (A34), (A35) and (A37). The results are shown in Fig. 2, for parameters appro- 
priate for the Crab pulsar: <r L = 3 x 10 4 , 7l = 200, (> y/aZ, so that the wind is initially supersonic) 
and lul/Q = 10 11 . The multiplicity parameter for these parameters is k = 8 x 10 3 . Solutions for 
three different initial conditions are shown, corresponding to starting points for reconnection at 
10, 10 3 and 10 5 x n,. Superposed on these solutions is the analytic asymptotic solution. For all 
values of the initial condition, the asymptotic solution is approached rapidly, and followed closely, 
until A ~ 1. The termination shock in the wind of the Crab pulsar is thought to lie at roughly 
r s ~ 10 9 tl. At this point, A ~ 0.01 and the energy flux is still dominated by Poynting flux. 



4.3. Validity of the solutions 

The range of validity of our treatment is limited by two factors. According to our solutions, 
the temperature in the current sheets decreases outward. The expression used for the gyro-radius 
in Eqs. (17) and (A32) assumes T > mc 2 , as does the choice 7 = 4/3 for the ratio of specific heats. 
This is not the case at very large radius and the range of validity is therefore restricted to 

7 < -g- (34) 

At the upper end of this range, A = 1/18, so the assumption A <C 1 is still valid. The corresponding 
radius T2 is 

TTUJ L 10 6 /X 30 . . 

r2 = 18^ rL = ^^ rL (35) 

At this point, the fraction of the Poynting flux transferred to the plasma is still small. For the 
Crab, one obtains (see Fig. 2) T2 ~ 10 10 rL, which already exceeds the radius of the standing shock. 
Thus, our treatment is valid up to the shock front, before which reconnection is indeed ineffective 
in the wind of this pulsar. 
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Fig. 2. — Results of a numerical integration the system of equations (A34), (A35) and (A37) for 
three different initial conditions, superposed on the analytic asymptotic solution (dashed line). 
The parameters used are those appropriate for the Crab: cjl = 3 x 10 4 , 7l = 200, ujl/Q = 10 11 . 
Plotted are the temperature T in the hot sheet, in units of mc 2 , the fraction a of particles in the 
sheet [An[ i /[An[ i + (1 — A)n(.]] the Lorentz factor of the flow 7 and the fraction of a wavelength 
occupied by the sheets A, as functions of the radius in units of the light cylinder radius. The initial 
conditions correspond to reconnection starting at r /Vl = 10 (solid line), 10 3 (dotted line) and 10 5 
(dashed-dotted line). In the Crab, the termination shock is located at log(r/ri,) ~ 9. 
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The second limitation concerns the assumption that the dissipation proceeds sufficiently quickly 
to keep the sheet thickness equal to the minimum value given by Eq. (12). This assumption holds 
provided the proper propagation time, r/(cj), exceeds the Larmor period. Because the sheet width 
is equal to the Larmor radius, this condition, 

- > 4vr 2 A7r L (36) 

7 

is equivalent, to within a factor of the order of unity, to the condition that pressure equilibrium 
within the sheet has time to be established. Substituting Eq. (31) and Z = 1/3, one finds 

7 < 7m = 3k. (37) 

If 7max > 7m, the flow accelerates only up to the point where 7 = 7m, after which no further 
dissipation can occur. This regime arises for 




- < A^ (38) 



2 . 10 3V^30 (3Q) 



or, alternatively, 



7max = ^7L > J ~^ (40) 

= 6-10 3 ^ (41) 

Otherwise, acceleration continues and the flow can, in principle, reach 7 max > unless it encounters a 
shock front. 

In the case of a pulsar moving through the interstellar medium, one can estimate the position 
of the termination shock by equating the magnetic pressure to the ram pressure of the medium: 

§ (W - ^ <«> 

where p is the density of the interstellar medium, and V the velocity of the pulsar through it. 
Taking V = 100kms~ and a particle density of lcm~ 3 , one obtains 

r s = 7xl0 14 ^cm . (43) 

Comparing this with the expression for r max , Eq. (14), one sees that the flow remains Poynting 
dominated up to the termination shock for all but the milli-second pulsars. 
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5. Conclusions 

We have examined the fate of a wave generated in a pulsar wind by the rotating, magnetized 
neutron star. This wave is built from the oscillating equatorial current sheet which, at large 
distances, may be considered locally as a sequence of spherical current sheets separated in radius 
by the distance nr^. As was suggested by Usov (1975) and Michel (1982), the wave decays 
because the particle number density eventually becomes insufficient to maintain the necessary 
current. Dissipation begins when the velocity of the current carriers reaches the speed of light or, 
almost equivalently, when the sheet width becomes equal to the Larmor radius (Coroniti 1990; 
Michel 1994). The dissipation process may be considered as reconnection of oppositely directed 
magnetic fields (Coroniti 1990). In the proper plasma frame, plasma from the inter sheet space 
slowly moves towards the sheet, which slowly expands, absorbing more plasma and magnetic energy. 
The distance at which the wave decays completely is proportional to the Lorentz-factor 7 of the 
flow [Eq.(9)], as found by Michel and Coroniti. In this formula, however, they inserted the initial 
Lorentz-factor, 7 ~ 100 — 1000 and concluded that in the case of the Crab pulsar the wave decays 
well before the wind reaches the termination shock. 

We find that the flow accelerates during the dissipation process. The reason is that in the 
freely expanding flow, hot plasma in the current sheet performs work on the flow. The restraining 
magnetic tension force is released by reconnection, but the accelerating pressure gradient remains. 
As a result, most of the magnetic energy is dissipated when the flow has accelerated to a Lorentz- 
factor which is of the same order as the maximal one. For the Crab, the corresponding distance is 
well beyond the standing shock, so we conclude that the wave does not dissipate before entering the 
shock. Using a simple estimate of the location of the termination shock for other isolated pulsars, 
we find this conclusion holds for all except those of milli-second period. In the Appendix, we show 
that our result applies not only to the singular current sheet structure discussed by Coroniti (1990), 
but also to a more general smooth distribution of current and magnetic field. 

Exactly how the wind energy dissipates remains a mystery, not only in the case of pulsars, but 
also in the closely related models proposed for gamma-ray bursts (Usov 1992, 1994; Blackman &: 
Yi 1998). At present, one can only speculate that dissipation might be possible in a combination of 
shocks and current sheets at the position of the bright equatorial X-ray torus observed in the Crab 
(Brinkmann, Aschenbach & Langmeier 1985). If the Crab is an oblique rather than a perpendicular 
rotator, a significant part of the energy flux at high latitudes is transferred by the axisymmetric 
part of the Poynting flux, which cannot dissipate by reconnection. It has been suggested that 
the release of this energy may be triggered by the kink instability (Lyubarsky 1992; Begelman 
1998), giving rise to the jet-like structure observed in the Crab Nebula, orientated, presumably, 
along the rotation axis of the pulsar (Hester et al. 1995). Thus, the idea that dissipation of the 
axisymmetric component and the wave component of the Poynting flux is fundamentally different 
has both observational and theoretical support. Current sheets are not expected to be responsible 
for the former. We have demonstrated in this paper that they also cannot be responsible for the 
latter, until the pulsar wind encounters the termination shock. 
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A. Equations of entropy-wave evolution in the two-timescale approximation 

A.l. MHD equations 

Consider a non-steady, axially symmetric, radial MHD wind. In spherical coordinates, the 
electromagnetic fields, fluid velocity and current density in the laboratory frame are B = B(r,t)(p, 
E = E(r, t)9, v = v(r, t)r, j = j(r, t)9, and the proper (i.e., in the fluid rest frame) energy density, 
pressure, temperature (in energy units) and number density are e'(r,t), p'(r,t) T'(r,i) and n'(r,t). 

Then the equations of MHD are that of continuity: 

J^nO + ^I^W) = (Al) 



(where 7 = 1/yl — v 2 /c 2 ), the energy equation (zeroth component of the divergence of the stress- 
energy tensor): 

|(^°°) + ^|:(^ 01 ) = (A2) 



where 



the entropy equation: 



T 00 = ( e '+p') 1 2 -p'+^-^- 

T 01 = (e' + p')7W^ (A3) 



1 (d P > VdrA r E _v B)j (A4) 



7 — 1 \ dt n' dt J V c 

(where d/dt = d/dt + vd/dr and 7 is the ratio of specific heats) and the two Maxwell equations 



(Ampere's law) and 



1 d , x IdE Air , , N 

- — {rB ) + -— + — j = (A5) 

r or c at c 



1 d , ^ N 1 dB . , . 

(rE) + -— = (A6) 

r or c ot 



(Faraday's law). The system is completed by Ohm's law 



j = o cl (E- V -B) (A7) 
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the ideal gas law 

p' = n'T' (A8) 

and an equation of state, for which we select 

p' = (7 - l)( e ' - n'mc 2 ) (A9) 

where m is the (mean) particle rest mass. In the following we will take 7 = 4/3, as appropriate for 
a relativistic gas. 

Eliminating n' from Eq. (A4) using the equation of continuity (Al) we find 

In the limit of a relativistic gas, e' — ► 3p', the set of equations (A2), (A10), (A5) and (A6) is then 
independent of n', which appears only in the continuity equation (Al). 



A. 2. Perturbative solution 

To solve these equations, we exploit the two-timescales present in the problem. The fast 
timescale is that of the pulsar rotation period if ast = P = 2ir/Q. We assume that the initial 
conditions in the wind close to the star have period P and that at any fixed radius, all quantities 
vary with this period. At a distance of a few light-cylinder radii (n, = c/Q), we assume that the 
wind has settled down into an almost stationary pattern in which the fluid speed does not vary 
on the fast timescale, although the density, pressure and especially the magnetic field do so. Our 
initial conditions at this distance, therefore, constitute an entropy wave comoving with the fluid. 
The wave evolves on the slow timescale which is the expansion timescale of the wind £ s iow ~ rP/rj^. 
For r » tl, we have conditions suitable for a two-timescale expansion, i.e., £f ast <C t s \ ow . The 
general procedure is: i) transform to fast and slow independent variables, ii) expand the dependent 
variables, iii) collect and solve the zeroth order equations, and iv) collect the first order equations 
and demand that the secular terms they contain vanish (Nayfeh 1973). 



First, we define a fast phase variable 



t 

'0 v * 



(All) 



where v w (r) is the speed of the pattern, which is to be determined. We now change variables from 
(r, t) to ((f), R), where the dimensionless 'slow' variable is defined as 

R = er/r L (A12) 

with £<1 and R ~ 1. In terms of the new variables, we find for the continuity equation (Al): 
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the energy equation (A2): 

8T 00 1 dT m e d 



d(f> v w d(f) R 2 dR 
the modified entropy equation (A10): 

, .d'y dp' , p' d . 

d(j) d(j) V w 0(f) 



+ ™^(-R i ^ 1 ) = (A") 



^f^>V) + 3^>-^ (A15) 



Ampere's equation (A5): 



Faraday's equation (A6): 



8E 1 dB e d ,__, 4vr . ,_„, 



dB 1 8E e d ,„„, ,. „, 

a* " ^ a* + fl«i (/;£) " ° ' (A1T> 

In these equations we have expressed all velocities in units of the speed of light (v — > cw etc.). 

We now expand the dependent variables, noting that for an entropy wave the zeroth order 
velocity is independent of the phase: 

v = v (R) + evi((f>, R) 
ri = n' (<p, R) + en[(<f), R) 
p' = p' o ( ( f ) ,R)+e P > 1 (0,R) (A18) 



etc. 

Substituting and collecting terms of order e° we find from the continuity equation (A13) that 
for a non-uniform wind (one in which 70^0 is a function of (f>) 

v w = v (R) . (A19) 

The Maxwell equations (A17) and (A16) then lead to 

dB _ 4tt^v 



n 



-k • (A20) 



From Faraday's law (A17), it then follows that the quantity Eq — vqBq is independent of (p. Fur- 
thermore, in order to describe a wave with both a reconnection zone and a region in which ideal 
MHD holds (a c — ► 00), it follows from Ohm's law Eq. (A7) that this quantity must be zero: 

E = v B . (A21) 
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The energy equation (A14) yields the pressure balance condition: 

2 




■OiU + A) - () ( A22 ) 



0<P V" U 8vr 7o 2 

Finally, using Eqs. (A19) and (A21) it can be seen that the entropy equation (A15) is satisfied 
identically to zeroth order. 

To first order in e, we have for the continuity, energy and Maxwell equations: 

Sl (nW - *£si*™*> <a23> 

d(j) \v L J d(p V w 4tt%v 

- WM^ V > < A24 > 

U E '~^) = -5^^>"^ < A25> 

I ( Bl " £*) = -535<***> ' <A26> 

In the entropy equation, one may substitute for the zeroth order current using equation (A20) to 
find 

470 d , 

1 (* " **)?£ = #£(*V*) + 3«| • (A27) 



47r 7 oV u ; d(j) R 2 dR y ,u u; ,u u 5i2 

The equations governing the evolution of the zeroth order quantities on the slow scale are given 
in the usual manner by imposing non-secular behavior on the first order equations (Nayfeh 1973). 
This ensures that the first order quantities do not grow to dominate the zeroth order terms of the 
expansion within 0{r/ri) wave periods. As frequently happens, the imposition of these regularity 
conditions suffices to determine the slow variation of the zeroth order quantities. Consider, for 
example, Eq. (A23), which is a linear, inhomogeneous equation for v\. The right-hand side is, by 
construction, a periodic function of <fi. Therefore, v\ will grow with (j) unless the integral of the 
right-hand side over a complete period vanishes. Applying these considerations also to Eq. (A24) 
leads to the conditions: 

AfjJ 2 7o „ /"clKJ = (A28) 



u^r^ 1 ) - &*** 



2tT / R 2 



#(eo+Po + 47r7 2 



o 
(A29) 
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Equation (A25) is needed only if the first order current is required, and Eq. (A26) yields only the 
conservation of the phase-integrated flux. Equation (A27) is integrated by parts to give, using 
(A26) 



2w Ad 

2w 



Q 



1 



47T70-R 



o 



27T Q 

dcpBo— (Rv B ) 



(A30) 



A. 3. Application to the striped wind 

For the striped wind (Fig. 1), the zeroth order solution is 



n' = n' h (R) 
P ' =p'(R) 
B = 
n' = n' c (R) 

B = B(R) 



> for < (p < ttA(R) and vr < (p < vr[l + A(R)] 



> for irA < 6 < 7r 



n' = n' c (R) 
B = -B(R) 



> for 7r[l + A(i?)] < 0<2tt 



together with the condition of pressure equilibrium between the hot and cold layers: 

B 2 



P 



87T7o 2 



B" 

"87 



(A31) 



Coroniti's estimate of the thickness of the neutral sheet gives 



A(R) 



P '{R)n 



n' h {R)eirB(R)v (R) 



and the ideal MHD condition outside the sheet implies flux freezing there: 



_d_ 
OR 



B(R) 



R l0 (R)n' c (R) 







(A32) 



(A33) 



This latter relation follows formally by using the first order form of the ideal MHD condition: 
E\ = voBi + v\Bq in Faraday's equation (A26) together with the continuity equation (A23) and 
the fact that in our present configuration Bq and no are independent of (f> outside of the sheet. 
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The five unknowns obey, in addition, equations (A28, A29, and A30). In the case of Eqs. (A28) 
and (A29) the integration over cj) may be performed immediately to give: 

^{R 2 l0 v [(l-A)n' c + An' h }} = (A34) 

^|fl 2 7o Vc 2 [(l-A)< + A<] +2 i? 2 7o 2 t, (l + A)|^| = (A35) 

[Eq. (20 and 22]. The left-hand side of Eq. (A30) is also straightforwardly integrated. On the 
right-hand side, however, we first rewrite the integration in terms of Bq\ 



4Art a ^.,.„, , ,.,.... a , A „ n _ i r,A^™i +Bl a {Rva) 



CJ 

(R 2 jov ) + 37o«o^ (Apo) 



R 2 dR \- wj ~i»»v dR \-ruj An 10 Rj 



2 OR v dR 



(A36) 



Now the integration over <f> can be performed unambiguously to give, using Eq. (A31): 

4AB' 2 d ,„, N „ d 



R 2 lovo)+ 3~f v — (AB' 2 ) + 



R 2 dR \- <u«uj ' "<»"" dR 

^K(.-^ + ^^4w^0. <A3 7 ) 

To integrate the entropy equation (A37) we first multiplying it by R?^qVq and combine terms 
with and without A to find 

AloVo M ( R2v ™ AB ' 2 ) ~ 2 ^ vR Jr ( Rv oloAB> 2 ) + — (R 2 v 2 ol 2 B' 2 ) = 0. (A38) 

The last term appears to be of zeroth order in A, but, because of the continuity equation (A34), 
it is in fact of first order. Substituting for B' using the flux freezing condition of Eq. (A33), this 
term can be reduced to 

2R 2 v olo n' c ^ R - (R 2 v 0l0 n' c ) = 2R 2 v 0l0 n' c —(R 2 v 0l0 A(n' c -n' h )) (A39) 

where the continuity equation (A34) was used. Returning to equation (A38), and using the substi- 
tution 

R 7 

v Ql0 A = — (A40) 

47TK 

where Z is defined in Eq. (18), one finds 

<\7 c\n' 

2Rn ' c dR + K(3Z " 1} + (3Z " 1)R IR = °" (A41) 

which integrates to 

(1 - 3Zf/ 3 R 3 n' c = constant . (A42) 



22 



The striped wind shown in Fig. 1 is a particular (and singular) idealization of a wind containing 
cold magnetized parts of opposite polarity separated by hot sheets. More generally, we can envisage 
an idealization in which n' , Bq and p' are all continuous functions of cf). First, assume Bq and n' 
take on the constant values B and n' c outside of the sheets. From the condition of pressure balance 

B 2 = Bl^ + ^lp'^) (A43) 

Defining the effective sheet width as 



B 2 Jo 2vr 
and the average particle density n' h within the sheet via: 



«-<)A = jT^ [n' o (0) - <] (A45) 

we find, from Eqs. (A28-A30) a system of equations which is identical to those obeyed in the 
singular idealization. 
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